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ABSTRACT 

Coherent synchrotron emission by particles moving along semi-infinite tracks is dis- 
cussed, with a specific application to radio emission from air showers induced by high- 
energy cosmic rays. It is shown that in general, radiation from a particle moving along 
a semi-infinite orbit consists of usual synchrotron emission and modified impulsive 
bremsstrahlung. The latter component is due to the instantaneous onset of the curved 
trajectory of the emitting particle at its creation. Inclusion of the bremsstrahlung 
leads to broadening of the radiation pattern and a slower decay of the spectrum at the 
cut-off frequency than the conventional synchrotron emission. Possible implications of 
these features for air shower radio emission are discussed. 

Key words: Plasmas-radiation mechanisms: nonthermal-cosmic rays 
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High energy cosmic rays can initiate an electromagnetic 
cascade in a medium where relativistic electrons and 
positrons can be produced in a volume with a longitudinal 
(along the line of sight) dimension being smaller than the 
relevant radio wavelength. So, particles form a coherent 
bunch, acting like a single charged particle that emits 
a short burst of coherent radio emission. Two radiation 
proces ses have been considere d: coherent Cerenkov radi- 
ation jArskar'van I Il962t ll965|) and coherent synchrotro n 
radiation dFalcke fc Gorham \\20oA. iHueee fc Falcke 1120031) . 
(Cohere nt radio emission from a ir sho wers was first c onsid- 
ered by iKahn fc Lerche I dl966ft . also IColgatd dl967D . but 
their theory was not explicitly based on geosynchrotron 
emission.) The former requires charge asymmetry, say an 
excess of electrons, and a relatively dense medium for 
coherent Cerenkov emission to be at radio frequency. For 
example, the Moon is a good target for high energy neutri- 
nos that can lead to a cascade in the lunar rocks. Excess 
electrons can develop leading to coherent Cerenkov emission 
at radio frequency with wavelength compara ble with or 
larger than the longitudinal size of the cascade dArskar'van I 
19621 . Il965"l. For air showers, coherent synchrotron emission 



is generally considered more i mportant than Cerenkov 
radiation (in the radio band) (Falcke fc Gorham I l2003t 



I Gorham. et al. 1 120041) . The shower produces a bunch of 
relativistic electrons and positrons emitting coherent syn- 
chrotron radiation in geomagnetic fields. Unlike coherent 
Cerenkov radiation, cohe rent synchrotron emission does not 
need charge asymmetry dFalcke fc Gorham 1120031) . 

So far, the relevant spectra of coherent synchrotron 
emission from air showers were commonly calculated 
using numerical sim ulation based on the re t arded- 
potential method llSuprun. Gorham fc Rosner I 120031 : 
IHueee fc Falcke 1 l2005al |b|). In this method , the radiatio n 
is calculated from the retarded potential \ Ja cksorJ Il99&f) . 
Although cohere nt synchrotron emission has been consid- 
ered analyticall^ j[Ajfjisif^fcBla. sT1l2002t IFalcke fc GorhanTI 
2003; Huc ge fc Falcke II2003T) . their calculation is based on 
the standard synchrotron radiation formula, which does not 
include the effect due to the particle's finite track. For air 
showers, effective coherent emission occurs in the core of the 
shower where most of radiating particles are created. Thus, 
it is of interest to consider radiation by a particle moving 
along a semi-infinite trajectory. Apart from the usual 
synchrotron emission there is emission due to the onset of 
the particle's curved trajectory. The latter component is 
referred to as the modified impulsive bremsstrahlung (MIB) 
as it is due to the combined effect of the usual impulsive 
(or prompt) bremsstrahlung due to particle (e ) creation, 
which is modeled as an abrupt jump in the particle's 
velocity from zero to c, and curvature of the particle's 
trajectory. It is worth noting that t he finite track effect was 
considered for Cerenkov radiation IITamm 1119391) and was 
taken into account explicitly in cal culation of cosmic-ray 
induc ed showers in a dense medium dZas. Halzen fc Stanevl 
11992ft . In the case of Cerenkov radiation, the finite track 
leads to a reduction in radiation intensity and modification 
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of angular distribution of the emission. Since the main 
objective of studying radio emission from air showers is 
to infer the properties of the cosmic rays that induce 
the showers, one needs to determine the radio spectrum 
accurately. In this paper, we present a general formalism 
for coherent synchrotron radiation from a nonstationary 
many-particle system, which takes account of MIB. The 
formalism developed here is based on t he single particle 
treatment iMelrose fc McPhedron"lll99ll) . in which radia- 
tion is due to a current associated with particle's motion 
in a medium. Here, the current is regarded as extraneous 
as it is different from that due to the plasma response to 
waves. The spectrum of radiation from a many-particle 
system is derived from the total current that is obtained by 
adding all the currents due to individual particles. In the 
relativistic limit as in the usual synchrotron radiation, the 
spectrum can be expressed in terms of the Airy functions 
and as a result, the radiation is highly beamed. 

In Sec. 2, a general formalism for synchrotron emission 
from a many-particle system is derived by including the ef- 
fect of MIB due to the effect of a particle's semi-infinite 
track. Coherent synchrotron emission is considered in Sec. 3 
and the application to air showers is discussed in Sec. 4. 



The orbit of a charged particle spiraling in a magnetic 
field B oriented along the z-axis can be written as 



x(t) = Xoc + -jy 

+v\\te x , 



sin (fit — tpo)e x + cos (fit — tpo)e y 



(3) 



where t 0, xo c is the initial position of the particle's gy- 
rocenter, tpo is the initial gyrophase, which is defined here 
as the the azimuthal angle of the particle's initial velocity 
relative to the magnetic field, fi = r/Q, e /'y, q = r\e, 77 is the 
charge sign, 7 = 1/(1 — u 2 /c 2 ) 1//2 is the Lorentz factor, and 
fi e = eB/m e is the gyrofrequency. The standard method 
to calculate the current is t o expand the exponential term 
in terms of Bessel functions ijMelrose fc McPhedro 'n"l ll99ll ). 
In the relativistic limit as in the case relevant here, instant 
emission can only be seen during a very short time interval 
At ~ Rg/cy 3 = /3x/fie7 2 , where R g = Di/fi is the gyrora- 
dius. Thus, the orbit can be expanded on tfi e ~ P±/"f 2 <C 1 
and the exponential in Q can be expressed into the form 

cot - k ■ x « d + ar + 6r 3 , (4) 

where r = t + to with to = fi -1 tan(<^ — tpo) and 



2 ENERGY SPECTRUM 

In cosmic ray air showers, electrons and positrons are created 
with a relativistic velocity, emitting synchrotron radiation in 
the geomagnetic field. Effective coherent emission by these 
secondary particles occurs at the core of the shower located 
very close to where most particles are produced. Thus, the 
finite track effect, in particular the initial position of the par- 
ticle's orbit, can be important and needs to be included ex- 
plicitly in the calculation of the synchrotron spectral power. 
In the following, we start with the single particle formalism. 



2.1 Single particle treatment 

Consider a charged (q) particle created at t = moving 
along a trajectory x(t) with a flight time T. The current 
associated with the particle is 



j(uj,k)—q / dt v(t) exp [i (cot — k ■ x(t))] . 



(1) 



where v — dx/dt is the particle's velocity. The current Q 
can be regarded as an extraneous current due to a single 
particle's motion (as compared to the induced current due 
to plasma response). In the usual application to radiation 
in astrophy sical plasmas, the time integration is taken from 
—00 to 00 ( Melrose 1986). A finite T, particularly the initial 
point at t = 0, introduces a boundary (finite track) effect 
into JTJ. The energy spectrum can be found from the ex- 
pression (iMelrose fc McPhedron"lll99lT) 



Um (k) 



1 , . ..a 

7T- \ e M ■ J\ , 

ISO 



(2) 



where e* M is the complex conjugate of the polarization ein 
of the wave emitted in the mode M. The single particle's 
spectral power can be derived from Pm = Um/T. 



tpo) + \ tan((£ — tpo) sin(0 — ^0) 



|fc±uj_fi 2 cos(0 - tpo), 



-k ■ xoc 



u — k\\v\\ 



2 f . v sin 2 (</>- tpo) 
3 cos((p — tpo) 



tan(0 — tpo). 



(5) 



We assume the observation direction is fc = k/k — 
(sin 8 cos (p, sin 8 sin <p, cos 8) and define spherical coordinates 
fc = e r . The projection of the current to the plane perpen- 
dicular to fc can be written as j ± = j — kk ■ j , which has 
the following components: 



erj 



v± 



s{4>-tpo) 



+i£lF' s'm((p — tpo) \ cos 8 — v^F s'md 

j± r p w ier/v^QF' cos(<p - tpo) e d , 
where the relevant integrals are defined as 



(6) 



F 



F' = 



dr e 



iiaT + bT-*) 



OF 
da 



AT T G ' . 



(7) 



The relativistic beaming (7 3> 1) implies that \<p — <C 1 
and \ol — 0\ <C 1, where a = arctan(wx/w|i) is the pitch angle. 
In these approximations, one has a w u)[l — n/3 cos(# — a)] ~ 
(o>/2)[2(l -n) + 7~ 2 + (8- a) 2 ] and b w (nwfi 2 /6) cos(0 - 
tpo) sin 8 sin a, where the refraction index is assumed to sat- 
isfy |1 —n\ < 1. 
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2.2 Spectrum 

One may write the energy spectrum as the energy radi- 
ated per unit frequency per unit solid angle, U(u>,k) = 
(oj 2 /8tt 3 c 3 ) J2 m U M (k), where U M (k) is given by © and 
the summation is made over two orthogonal modes. One 
finds 



U(u>,k) = 



167T 3 £oC 3 



\j±e\ +\J±<i>\ 



167T 3 £oC 



(a- 9y\Ff + -f sin 



(8) 



The integrals in (0 can be expressed in terms of the Airy 
functions provided that the flight time T is much longer 
than the duration of the synchrotron pulse (At ~ l/f2 e 7 2 ). 
In this approximation, one may take the limit T + to — * oo 
and the integrals reduce to the form 



F' 



(3&) 1 / 3 



(36) 2 / 3 



Ai(0+iGi(0 



Ai' (0+iGi'«) 



d${a,b) 
da 



lj \ 2 / 3 2(I-7i)+7- 2 + (6»-a) 2 



(3&)V3 



(2^) 



i/)o) sin sin a] 1 / 3 



dr e 



i(aT + 6r ) 



(9) 



(10) 



(11) 



(12) 



where Gi«) = (l/3)Bi(0-£[Ai(OBi(0- Ai(£)Bi(0]<, 
Ai(£) and Bi(£) are the Airy functions (see Abramowitz & 
Stegun 1970). Similar to Ai(£) and -Ai'(f), both Gi(£) and 
— Gi'(£) are a decaying function for £ ^ (as shown in 
figure . The first terms in both square brackets on the 
right-hand sides of @ and 110H correspond to synchrotron 
radiation for a particle moving along a semi-infinite trajec- 
tory. The terms Gi and Gi' describe MIB arising from the 
combined effect of the semi-infinite track's initial point and 
curvature. One can show that MIB has features of usual im- 
pulsive bremsstrahlung (cf. Eq |161 . i.e. emission due to an in- 
stantaneous change in the velocity from zero to v ~ c at t = 
(cf. Eq. 1) (lLandau fc Lifshitz ll97ltlGrichine 12001 . One 
should emphasize here that the impulsive bremsstrahlung 
considered here is different from the usual bremsstrahlung 
by a charged particle interacting with the Coulomb field of 
nuclei in matter. It can be shown that the last terms ($ 
and dQ/da) on the right-hand side of both J§J and HOI 
can be ignored provided that \<p — tpolf <C Tr(2Q, e 'y 2 /uj) 1 ' 3 . 
The characteristic frequency can be estimated from a 3 ~ b, 
which leads to u ~ 7 2 f2 e . 

The spectrum can be written as sum of usual syn- 
chrotron emission (Usyn) from a semi-finite track and MIB 
emission (Ub) due to the velocity jump at t = and the 
trajectory's curvature, that is, 

U = U syn + U h , (13) 
with 
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Figure 1. Plots of Gi(£) (solid) and Gi'(£) (dashed). 
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f2 e sin a 
ay 



(*£ a Ai'(£) 



(14) 



U h 



fi e sin a 

a~f 



(VGi'(£) 



(15) 



where we neglect both the term (1121 and its derivative. 
It is often convenient to separate the radiation into two 
components, with one polarized perpendicular to the orbit 
plane, corresponding to the first terms in 1141 and 1151 . and 
the other polarized in the plane, corresponding to the sec- 
ond terms involving derivatives. Specifically, one may write 
U syn = C4yn_L + Usyn|| and Ub = Ub± + U b \\- Notice that 
Usyn is smaller by a factor of 4 than that for normal syn- 
chrotron emission from a full infinite track (cf . Eq. 1171 . Tak- 
ing the limit £l e — > one can easily verify that Ub reduces to 
the familar form for the prompt bremsstrahlung in the case 
where a particle's velocity abruptly changes from zero to a 
constant velocity jLandau fc Lifshitz Ill97ll) . In this limit, 
the angle a is re-interpreted as an angle that the velocity 
makes with respect to the z axis. Since one has £ — > oo and 
hence 7r£Gi(£) — » 1 (Abramowitz & Stegun 1970), one finds 



U b (0J,k) : 



47r 3 £oC 



7-2 + 1 



(16) 



which is the same as that for the case where a particle in- 
stane ously acquires a constant velocity llLandau fc Lif shitz 
11971ft . 

Apart from its property of the usual impulsive 
bremsstrahlung, MIB also has a feature of synchrotron ra- 
diation (cf. Sec. 2.3), i.e. Ub increases with frequency as a 
power-law, similar to the usual synchrotron spectrum at low 
frequencies. Such similarity can be understood as due to that 
the emitting particle's trajectory is curved, while for the 
usual impulsive bresstrahlung the particle's trajectory is a 
straight line. Because of this similar feature to synchrotron 
radiation MIB is nonzero when one applies it to pair cre- 
ation. 
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Figure 2. Plot Uu (arbitrary scale) as a function of (9 — a)j 
with 4> = ipo f° r a single particle. The dashed line corresponds to 
the usual synchrotron emission with Un = 4[/ syn || , in the absence 
of MIB emission. We assume 7 = 80, n = 1, the gyrofrequcncy 
!! e = 5x 10 6 b _1 , u>/2tt = 100 MHz and the pitch angle a = tt/4. 
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Figure 3. Plot U± (arbitrary scale) as a function of (9 — a)j 
with <j> = ?/>o- The parameters are as in figure [?""| Apart from a 
reduction in intensity, the modified synchrotron emission (usual 
synchrotron plus MIB) has a much wider profile than the usual 
synchrotron emission. 



2.3 Usual synchrotron radiation 

The conventional synchrotron emission can be reproduced 
by adding the other semi-infinite orbit from —00 to in 
Specifically, one first drops the terms $ in @ and 
<9<£>/<9a in 110II . both of which are cancelled out by the 
corresponding terms from the negative semi-infinite trajec- 
tory, and then replaces F and F' in J""} with 2Ke(F) and 
2Re(F'), respectively. Thus, the bremsstrahlung terms can- 
cel o ut. One then reproduces the usual synchrotron for- 
mula ijMelrose fc McPhedron"lll99ll) 



U(uj,k) 



4t/ s 



syn 

2 



V Q, sin a J 



7 - 2 + {9 -a] 



7 + 1 



'-a)' 



:K\ /3 {p) +Kl /3 {p) 



(17) 



where n — 1 and p = (2/3)£ 3//2 . We rewrite the Airy func- 
tions in terms of the modified Bessel functions and the ap- 
proximations \Q — a\ <C 1 and \<j> — tpo\ <C 1 are used. 

The angular distribution of single particle's spectrum 
1131 is shown in figures 12*151 Here, the radiation is separated 
into the parallel component, Uu = £4^11 + Uwi and perpen- 
dicular component, U± = U syn ± + Uh±- The inclusion of the 
boundary effect leads to an overall reduction in intensity and 
a significant broadening of the angular profile. Figures[f""|and 
[7|show a comparison of the two components Ub and f/ syn as 
a function of frequency. These two components are similar to 
each other at low frequencies. However, !7b|| has two cutoffs, 
with the lower one determined from the zero of Gi'(£) (cf. 
figure ["[J. For the perpendicular polarization, Ub± levels out 
at high frequencies and behaves much like the usual impul- 
sive bremsstrahlung. The energy spectrum for the parallel 
polarization is shown in figure |H] The synchrotron emission 
drops off exponentially above the critical frequency. Since 
Gi(£) ~ l/£ and Gi'(£) ~ 1/C 2 > which drop off much slower 
than the exponential decay of Ai(£) and Ai'(£) at a large 
£, the bremsstrahlung component decays much slower than 
the usual synchrotron emission. 
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Figure 4. Angular profile (arbitrary scale) as in figure [2] The 
solid and dashed lines correspond respectively to U^n and (7 syn ||. 
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Figure 5. Angular profile (arbitrary scale) as in figure El The 
solid and dashed lines correspond respectively to Ub± and U ayn ±. 
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Figure 6. Energy spectrum (in JHz^sr" 1 ) as a function of 
frequency u>/2tt (Hz) for </> = tpo for a single particle. The solid 
and dashed lines represent 27rt/i,|| and 27rC/ syn ||, respectively. We 
assume (6-0)7 = 0.5, N = 10 7 , 7 = 80, n = 1, Q e = 5x 10 6 s" 1 , 
(8 — 0)7 = 0.5, and tp = ipo- At low frequencies, the spectrum is 
very similar to that for the usual synchrotron, i.e. U <x a; 2 / 3 for 
u) <C u c ~ fi e 7 2 . 



Figure 8. Energy spectrum 2-nUn 



2n(U h 



|) (m 



J Hz - 1 sr — 1 ). The dashed line corresponds to synchrotron emis- 
sion 87r[/ syn ||. The parameters are as in figure IB1 In comparison 
with the normal synchrotron emission (dashed line), due to the 
MIB the spectrum (solid line) drops off much slower than the 
exponential decay at u> 2> u) c . 
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where N is the total number of charged particles. The co- 



10 10 10 

m/2% (Hz) 

Figure 7. Energy spectrum as in figure [S] for the _L-polarizcd 
components. Notice that the spectral component 2-kV\ 3 x_ (solid 
line) tends to level out at high frequencies where the approximate 
form 1161 applies. 



2.4 Many-particle system 

The single particle formula can be extended to a many- 
particle system by adding all the currents from individ- 
ual particles. Let the sth particle be created at time t s . A 
nonzero t 3 adds a phase term u>t 3 to ©. The total current 
can be obtained summing up all individual currents given 
by JSJl with all relevant quantities labelled by s. Then, the 
total energy spectrum is derived as 



[/tot (u, k) 



167T 3 £oC ^—^ 



-I sin a s sin a 3 i r s t s i 

Is 7s' 



r? s rv (,a s - e)[a s , - e)F s F* s , 
'»', (18) 



herence is determined by the phase (p ss r = d s 
by 



given 



tp ss > w (j(t 3 — t s i) - k • (x 0s — x 0b i) 

-— [1 - nf3 s cos(6> - a s )] {</> - ipos) 

1 l s 

+■£- [1 - np s , cos{6 - a 3 ,)] (<f> ~ Voa')- (19) 

We use the following expression for the initial position: 
xo = xqc — («x/f2)(sin tpo e x — cos ipoe y ). The usual sponta- 
neous emission corresponds to that the phase is \tp S s'\ ^ 1 
for s 7^ s' and that only terms of s = s' contribute to the 
total spectrum. Therefore in the case of spontaneous syn- 
chrotron radiation, the total spectrum can be written as 
CAot = NU, where U is the single particle's spectrum aver- 
aged over the particles' momentum distribution. Since the 
initial gyrophase does not enter the final form of the energy 
spectrum, spontaneous synchrotron radiation is axially sym- 
metric with respect to the magnetic field line direction, i.e. 
the angular pattern of emission depends only on 6 not (j>. In 
the case of coherent synchrotron emission (cf. Sec. 3), such 
symmetry is broken since the maximum coherence depends 
explicitly on the initial gyrophases (cf. Eq. I19H . 



3 COHERENT SYNCHROTRON EMISSION 

Coherent emission occurs provided that the majority of 
emitting particles satisfy the condition \<p s — tp s r\ <C 1. In 
general, one can calculate the total spectrum numerically 
using 1181 for a given distribution of the particle injection 
time, position and initial velocity. In some special cases, one 
may write down its analytical form. The simplest case is 
where all particles have the same initial velocity, which is 
considered here. 

One may write the total spectrum into the form 
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Uu 



N 



S±(u)U ± + S {[ (ui)U ll 



(20) 



where 1 ^ Su t ±(ui) ^ N is called the coherence factor. 
The value S\\ t ±(cu) = N corresponds to completely coherent 
emission and S^±(lu) = 1 to spontaneous emission. When 
the number densities of electrons and positrons are equal, 
j±e does not contribute to the coherent power as contribu- 
tions from electrons and positrons cancel out. This leads to 
S± = 1. The total spectrum is 



Utat 



N U±+S\\(u)U n 



167TEocsin a 



V Q. e sin a / 



2/3 



Ai' 2 (C)+Gi' 2 (0 



(21) 



(22) 



where £ « (aj/2r2sina) 2/3 [7^ 2 + (0 - a) 2 ] and n = 1 is 
assumed. In contrast to coherent Cerenkov emission, which 
requires a charge asymmetry (net charge) (Askar'yan 1962, 
1965), coherent synchrotron emission can occur for a neutral 
plasma. In the case of charge symmetry, the polarization is 
linear, perpendicular to the plane of the magnetic field and 
wave vector. The coherence factor is then given by 



i 

iV 



N 

E 



COS [ul(t s — t s r) — k ■ (xqs — Xqs')] ■ (23) 



Eq. IlL'Mli is sum of phasors, which can be modeled as (e.g. 
Hartemann 2000) 

w 1 - (cosO) 2 + 7V<cose) 2 

r i V 2 

+2(cos9) N(l - (cos9) 2 ) , (24) 



where N 3> 1, B = ui(t s — n ■ xos/c), n 
average is made over a distribution P(e), 

(cos e) = / de p(e) cos e. 



kc/cu, and the 



(25) 



For a uniform distribution P(Q) = 1, corresponding to spon- 
taneous emission, one has S\\(u>) = 1, and for P(<d) = 5(0), 
one has Sm (ui) — N, corresponding to completely coherent 
emission. As an example, we consider the case in which par- 
ticles are injected at the same time, say at t s = 0, with 
a gaussian profile in the longitudinal (along fe) spatial dis- 
tribution with a width Al. The probability can be written 
as 



P(G) 



-(e/Aey 



0FAe 



(26) 



where Ae = uiAl/c. Then, one finds (cose) = 
exp(-A6 2 /2) = exp[-(cjAZ/c) 2 /2]. A plot of S\\ as a func- 
tion of Al is shown in figure [5] Similarly, if all particles are 
injected at the same location with a gaussian profile with a 
width At in time, one has (cose) = exp[— (ojAt) 2 /2]. This 
shows that to attain effective coherence, the spread in time 
of particle injection must be At ~ The spectrum 

from a gaussian bunch is calculated using 121 H and is shown 
in figure [TO! The spectrum consists of two regions, separated 
by a critical frequency u> co h ~ (21n7V) 1/2 c/AZ: the coherent 
(Su > 1) emission region u < u) co h, where the polarization 
is predominantly linear, and spontaneous (S\\ = 1) emission 
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Figure 9. Coherence factor vs bunch size. The solid and dashed 
lines correspond respectively to a gaussian bunch (1261 and a 
bunch of T-form 1271 with A = 5. One assumes the total number 
of emitting particles TV = 10 7 . The bunch size is in units of the 
wavelength A. 



region u) ^ (jj co h, where the polarization is elliptical as both 
Un and U± are present. The calculation of the coherence 
factor can easily be extended to other types of distribution. 
In particular, the F-probability dist ribution is thought t o 
be the more relevant for air showers (lAgnetta et al. Ill997f) . 
which is discussed in detail in Sec. 4. 

As the secondary particles from air showers have a dis- 
tribution in momentum, one needs to use the full expression 
1181 1 to calculate the spectrum. When the particle's momen- 
tum distribution is included, an analytical form for the spec- 
trum can be derived only in some special cases. For example, 
when the bunch size is smaller than the relevant wavelength, 
one may ignore the phase term 1191 and expresses the spec- 
trum as a integration of the single particle formula over the 
particle's distribution. The total intensity is about JV 2 times 
larger than that by single particle. 



4 APPLICATION TO AIR SHOWERS 

The formalism developed in the previous sections can 
be applied to radio emission from extensive air show- 
ers (EAS). Detailed modeling requires numerical model- 
ing of EAS, which was already consider e d by several au- 
thors llSuprun. Gorham fe Rosner 1 120031 : iHueee fe Falckel 
l2005al lbl). Here we only discuss qualitatively implications of 
the modified synchrotron emission for radio emission from 
EAS and further detailed modeling will be considered else- 
where. A notable feature of the modified synchrotron emis- 
sion is that the corresponding intensity derived from a semi- 
inifinite orbit is smaller than the conventional synchrotron 
emission by about a factor of 4. In principle, this feature 
can be tested against observations provided accurate radio 
spectra with well calibrated intensities become available, 
for ex ample with the planned Low Frequency Array (LO- 
FAR) jFalcke fe Gorham lEofB) . Other features include an 
angular broadening of the radiation pattern and a slower de- 
cay of the spectrum above the critical frequency (to ~ Q e j 2 ). 
The broadening due to prompt bremstrahlung is more pro- 
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(B/2it (Hz) 

Figure 10. Spectrum 2nUtot (in J Hz - 1 sr — 1 ) from a gaus- 
sian bunch (solid line) with a longitudinal size AZ = 1 m and 
the T-pdf bunch (long-dashed line) with A = 3. The short- 
dashed line corresponds to the spectrum derived by ignoring 
the prompt bremstrahlung. The parameters are as in figure [B] 
There is a critical frequency cj co h that separates the spectrum 
into the coherent region ui < £J GO h and spontaneous emission re- 
gion ui ^ i^coh- For the gaussian bunch, the critical frequency 
is u> coh /2iT ~ (21nAf) 1 / 2 c/2vrA/ ~ 270 MHz and for the T-pdf 
bunch it is ui coh /2n ~ {c/ Al)N 1 / < ^+ A ') m 1 GHz. 



nounced for the perpendicular polarization than for the par- 
allel polarization. Since the perpendicularly polarized com- 
ponent depends on the charge excess, such broadening is 
important when there is significant charge asymmetry. The 
spectral hardening occurs near the synchrotron cut-off fre- 
quency, which is much higher than the transition frequency 
LUcoh- So, this spectral feature may not be observable for air 
showers. 

In the following we derive the coherence factor using 
the procedure in Sec. 3. As an example, one assumes that 
an air shower occurs at a few km height and that all emit- 
ting particles are located in the shower maximum. As the 
zeroth order approximation, the near field effect can be ig- 
nored and the energy spectrum (Eq. 11811 derived in Sec. 2 is 
applicable. The near field effect may need to be included if 
the maximum of an air shower develops very near the ground 
level. The re are extensiv e discussions of air showers in the 
literature (Gaisser 1990) and in principle, one can obtain a 
distribution of particle's injection time (t s ), initial momen- 
tum (yvo,,) and position (airjs)- Assuming the primary cos- 
mic ray energy to be E p , the number of secondary electrons 
and positrons can be estimated as N ~ E p / , ym e c 2 ~ 10 7 
for E p ~ 10 15 eV and 7 = 80. In the practical situation 
these particles are injected over an extended range rather 
than at a single fixed height. For EAS, the distribution in t s 
for secondary electrons/positrons can be inferred from mea- 
surements of arrival times of charged particles (muons plus 
electrons/positrons). One should note that these measured 
times are not actual arrival times of electrons and positrons 
since muons arrive earlier than electrons/positrons. The ar- 
rival times can b e modeled by a F-probability distributio n 
function (T-pdf) jAgnetta et al. ll997l : fAntoni. et al. l200ll) 
Thus , the probability can be written as (iHueg^^^ickel 
2003): 



\ (*X 

Aer(i + A) \aq) 



-8/A9 



(27) 



for 9 > and P(6) = for 6 < 0, where F(x) is the Gamma 
function and the power index A can be estimated from the 
distribution of electrons arrival times. The standard devia- 
tion is cr r = (1 + A) 1/2 AB. From (2D, 

one obtains 



(cose) = (1 + Ae 2 

(0) = (1 + A)A0. 



-(l+A)/2 



{1 + A) arctan A9 



(28) 
(29) 



The corresponding critical frequency is given by w co h ~ 
(c/AtyN 1 ' < - 1+A '. This frequency defines a transition from 
coherent to incoherent emission. For AZ = 1 m, N = 10 7 
and A = 5, one has u) co } x /2tt ~ 1.2 GHz. The coherence fac- 
tor can be obtained by substituting 1281 1 for 1241 1. Here we 
write down the two limiting cases: 



N(l + Ae 2 



(1 + A) arctan A6 



w < c^coh, (30) 



UJ 3> U>coh- 



The dashed line in figure [5] shows S\\ as a function of 
the bunch size AZ, featuring a higher transition frequency 
w co h than the gaussian bunch. The low-frequency approx- 
im ation in <|30[l correspond s to the nonsquared form given 
bv iHueee fc Falcke I ll2003l) except that the cosine factor is 
retained here. 

The coherent emission can be elliptically polarized if 
there is a charge excess as it is the case for air showers. 
Assuming the excess number of electrons is N c , one has 



S± 
Sn 



Nc 
N 



(31) 



The right-hand side is sensitive to the cut-off energy of the 
excess electron energy. For a cut-off near the MeV energy, 
the excess of electrons can reach N c /N ~ 10 — 15%. (Note 
that cascades in a dense medium such as rocks can give rise 
to about N c /N ~ 20% excess negative charge.) As a result, 
the coherent emission can be elliptically polarized with an 
ellipticity ~ (5j_/S||) 1/2 (5c = {N C /N) 2 S C ~ 0.38 c , where S c is 
the typical ellipticity of single particle's radiation. 



5 CONCLUSIONS 

Synchrotron emission by a particle moving along a semi- 
infinite trajectory is considered. Since effective coherent syn- 
chrotron emission by secondary particles in an air shower 
occurs in the core of the shower where most emitting parti- 
cles are created, the initial point of the particle's trajectory 
need be included explicitly. It is shown that radiation from a 
particle moving along a semi-infinite track can be separated 
into the usual synchrotron emission and the bremsstrahlung- 
like emission (MIB). The latter is due to emission as the 
result of onset of the particle's curved trajectory. The spec- 
tral intensity of the modified synchrotron emission (usual 
synchrotron plus MIB) is lower than that for normal syn- 
chrotron emission, roughly by a factor of 4. It is interest- 
ing to note that such reduction is consis tent with the re- 
cent result from numerical simulation by iHueee fc Falcke 1 
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(2005b). The radiation pattern has a broader angular dis- 
tribution than the usual synchrotron emission. This feature 
is especially pronounced for the perpendicularly polarized 
component. The spectrum has a much slower decay above 
the critical frequency uj ~ 7 2 f2 e , while the usual synchrotron 
spectrum has an exponential cutoff above the critical fre- 
quency. In the application to radio emission from air show- 
ers, the reduced intensity can be verified in principle pro- 
vided accurate observations of the radio spectrum are avail- 
able. Although there exist early observations of air shower 
radio emission, there are some uncertainties in determinin g 
the calibration factor iAllanll97l] : lAtrashkevich et allll978fi . 
The current LOFAR Prototype Station (LOPES) and the 
future LOFAR may provide a better opportunity to test the 
predicted spectrum. Since the transition frequency (to spon- 
taneous emission) is much lower than the cut-off frequency, 
change near the cut-off frequency may not be observable. 
The broadening of the radiation pattern occurs mainly for 
the perpendicularly polarized component, which may be ob- 
servable provided that there is significant charge asymmetry 
in the emitting plasma. 

A major advantage of the formalism presented here 
is that the initial conditions including the time of parti- 
cle creation, initial velocity and gyrophase all appear in 
the phase 11911 and these quantities can be modeled sta- 
tis tically. Although the sin gle-particle formalism was used 
bv lWhite fc Melrose I (I1982T) to treat coherent gyromagnetic 
emission, in their calculation, the finite track effect was not 
considered. It is shown here that the distribution of the par- 
ticle injection time is important in determining the coher- 
ence. The semi-infinite track approximation adopted here is 
valid provided that the emitting particles are highly rela- 
tivistic with 7 3> 1. In the relativistic limit, the synchrotron 
pulse duration (l/fi e 7 2 ~ 10~ 8 s) is much shorter than the 
typical flight time T ~ 10~ 6 s for a particle's free path 500 
m, and therefore the orbit can be approximately regarded as 
semi-infinite. If electrons and positrons from an air shower 
have an energy cutoff extending to MeV energies with 7 ~ 1, 
the synchrotron approximation is no longer valid and one has 
cyclotron emission instead. In this finite orbit needs 

to be considered. One may extend the calculation in Sec. 3 
and 4 to include the particle distribution in momentum and 
this requires a numerical approach, which is not considered 
here. 
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